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FOREWORD 


This  report  summarizes  the  results  of  a  study 
undertaken  by  Dr.  Bernard  Budiansky  for  LMSD 
Research,  Lockheed  Aircraft  Corporation.  The 
vork  was  performed  tinder  Air  Force  Contract 
No.  AF-04(6U7)«l8l.  Dr.  Budiansky,  consultant 
to  Lockheed  during  the  summer  of  1953,  is 
Associate  Professor  of  Structural  Mechanics  at 
Harvard  University. 
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ABSTRACT 


Theoretical  calculations  art  made  of  the  natural  modes  and  frequen¬ 
cies  of  small-amplitude  sloshing  of  liquids  in  partially  filled 
circular  canals  and  spherical  tanks.  An  integral-equation  approach 
is  used  to  analyze  the  circular  canal  for  arbitrary  depth  of  liquid. 
A  similar  approach  for  the  spherical  tank  provides  modes  and  fre¬ 
quencies  for  the  nearly- full  and  half- full  cases.  These  results, 
together  with  the  known  behavior  of  the  nearly-empty  tank,  are  used 
in  conjunction  with  the  trends  established  for  the  circular  canal 
as  a  basis  for  estimating  frequencies  for  arbitrary  depth  of  liquid 
in  the  spherical  tank.  The  dynamic  analysis  of  the  container-fluid 
system  by  means  of  the  mode-superposition  approach  is  discussed,  and 
modal  parameters  required  in  such  analyses  are  presented. 
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NOTATION 


geometrical  coordinates  (see  Fig.  2-1) 
nondimensional  coordinates 
canal  or  sphere  radius 

geometrical  parameters  related  to  liquid  depth 

velocity  distributions  in  x,  y;  and  z  directions,  respectively 

oscillation  frequency 

time 

velocity  potentials 
displacement  potential 

acceleration  of  free  fall  relative  to  container  vail 

frequency  parameter;  \  =  cs^R/g 

liquid  depth  (variable) 

complex  variables 

complex  potential 

kernel  function  (circular  canal) 

kernel  function  (spherical  tank) 

mode  shape  for  circular  canal 

mode  shape  for  spherical  tank 

^  P  *(p) 

cylindrical  coordinate  system 
nondimensional  coordinates 
p/P 
PP 

complete  elliptic  integrals  of  first  and  second  kind,  respectively 
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modal  parameters 
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mass  of  container 
density  of  liquid 
mass  of  liquid 

external  transverse  force  acting  on  container 
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3ec*?cn  1 
INTRODUCTION 

For  sufficiently  small  amplitudes  of  motion,  the  dynamic  effects  of  the 
sloshing  of  a  liquid  in  a  partially  filled  container  can  be  analyzed  in  terms 
of  the  natural  inodes  and  frequencies  of  the  small,  free-s.,rface  oscillations 
of  the  liquid.  This  paper  is  concerned  with  the  calculation  of  the  modes 
and  frequencies  of  partially  filled  circular  canals  and  spherical  tanks, 
with  attention  restricted  to  the  modes  that  would  be  induced  by  lateral 
accelerations  of  the  containers.  The  usual  assumptions  of  inviscid,  in¬ 
compressible  flow  are  made,  and  the  container  is  considered  to  be  rigid. 

The  dynamic  analysis  of  the  container- fluid  system  by  means  of  the  mode- 
superposition  approach  is  discussed,  and  modal  parameters  required  in 
such  analyses  are  presented. 

The  case  of  the  half- full  circular  canal  was  considered  by  Rayleigh  (see 
Ref.  1,  p.  444),  who  carried  out  an  approximate  energy  solution  for  the 
fundamental  frequency.  Graham  (Ref.  2)  estimated  the  fundamental  frequency 
of  the  half-full  sphere  by  conceptually  dividing  the  liquid  into  a  collection 
of  half-disko  and  then  using  Rayleigh's  result  in  an  unspecified  averaging 
process . 

The  frequency  and  mode  calculations  of  this  paper  are  based  on  a  rigorous 
integral  equation  approach.  In  the  case  of  the  circular  canal,  the  appro¬ 
priate  integral  equation  is  fonnulated  explicitly  for  an  arbitrary  depth  of 
liquid,  and  numerical  calculations  of  the  first  three  pertinent  modes  and 
frequencies  are  performed  for  a  range  of  depths  by  means  of  a  matrix 


1-1 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILE  SYSTEMS  DIVISION 


LMSD-5151 


approximation  to  the  integral  equation.  For  the  spherical  tank,  rigorous 
analyses  are  made  for  only  the  nearly-enrpty,  half-full,  and  nearly- full 
tanks;  results  for  the  intermediate  depths  are  estimated  by  interpolation 
and  by  means  of  the  trends  established  for  the  circular  canal. 

The  results  for  the  fundamental  frequency  for  the  spherical  tank  are  compared 
with  available  experimental  data. 
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Section  2 

SLOSHING  MODES  AND  PREQUENCIES 

2.1  BASIC  DIFFERENTIAL  EQUATIONS 

The  differential  equation  and  boundary  conditions  governing  free- surface 
oscillations  (Ref.  1,  p.  3 63)  will  be  displayed. 

The  liquid  is  supposed  to  be  bounded  by  a  horizontal  free  surface  F  and  a 
rigid  wetted  surfe.ce  S  (see  Fig.  2-1),  and  small,  hannonic  oscillations 
are  assumed.  The  velocities  in  the  x,y,z  directions,  respectively,  are  then 


u  (x>  y>  z)  Bin  cut. 
v  (x,  y,  2)  sin  cob. 
x  (x,  y,  z)  sin  tut. 


(2.1) 


where  ©  is  the  circular  frequency. 


Since  the  flow  is  irrotational,  there  exists  a  velocity  potential 


♦  (x*  y>  z)  3in  cut. 


such  that 


(2.2) 
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The  continuity  ecjuation  for  incompressible  flow  is 

V2*  «  0 


(2.3) 


The  condition  of  tangential  flow  at  the  rigid  surface  requires  that 


b* 


»  n  =  =  0  on  8 


(2.4) 


where  n  is  the  unit  normal  vector  to  S. 

Finally,  the  condition  of  constant  pressure  at  the  free  surface  together 
with  the  linearized  Bernoulli  equation  leads  to  the  condition 


bt  (  m2 


♦  on  F 


(2.5) 


where  g  is  the  acceleration  of  free  fall  relative  to  the  container. 


Equations  (2.>),  (2.4),  and  (2.5)  govern  the  eigenvalue  problem  for  the 
frequency  a>.  Note  that  (2.1),  (2.2),  and  (2.5)  imply  that  the  shape  of  the 
vertically  displaced  free  surface  during  oscillation  is  given  by  ♦(*,  0,  z). 


The  analogous  two-dimensional  problem  of  transverse  sloshing  in  a  uniform 
channel  is  formulated  in  an  obvious  fashion  by  suppressing  the  z  coordinate. 

2.2  CIRCULAR  CANAL 


The  geometrical  parameters  for  the  partially  filled  circular  canal  are  shown 
in  Fig.  2-2}  the  depth  of  liquid  is  measured  by  e  .  which  varies  from  -1  for 
the  empty  canal  to  1  for  the  full  vessel. 

It  is  clear  that,  because  of  geometrical  symmetry,  the  natural  modes  involve 
free- surface  vertical  displacements  that  are  either  symmetrical  or  antisymmet- 
rical  about  the  orgin.  Since  the  symmetrical  modes  would  not  he  induced  by 
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transverse  motion  of  the  canal,  attention  is  restricted  to  the  antisymmet¬ 
ries!  modes. 


The  procedure  for  setting  tip  an  integral  equation  governing  the  antisynaet- 
rical  modes  follows. 


Consider  the  fluid-filled  container  consisting  of  the  original  wetted  boundary 
P-lus  its  reflection  about  the  |  axis  (Pig.  2-4).  Place  a  two-dimensional 
sink  of  strength  2w(j)d5  at  £  and  a  source  of  equal  strength  at  ,  and 
smear  out  both  over  a  strip  of  length  d£  .  Denote  the  resulting  velocity 
potential  for  the  internal  flow  by 

v(£)  •  U,  tjj  I)  dl 

Since  the  velocity  potential  is;  arbitrary  to  within  a  constant,  it  will  be 
assumed  that  ®(C,  tjj  |)  =  C,  and  hanse. *»  -  $(“§>  I)*  By 
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Fig.  2 


Fig.  2- 
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3  Ifondimensional  Geometrical  Parametero 
and  Coordinates  for  Circular  Canal 
and  -Spherical  Tank 


4  Location  of  Source-Sink  Pair  in 
Symmetrized  Domain 
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symmetry,  the  vertical  velocities  along  r{  =  0  vanish  everywhere  except 
on  the  strips  of  length  dg  at  5  and  -|  ,  where  just  belov  the  (  axis 
they  are  v(g)  and  -w(g)  ,  respectively.  It  is  now  possible  to  construct 
the  velocity  potential  for  the  free- surface  oscillation  of  the  configuration 
of  Pig.  2-3  by  means  of  a  distribution  of  sources  and  sinks  along  the  {  axis* 
It  is  necessary  only  to  satisfy  the  free-surface  condition  (Eq.  [2*7]) 
relating  the  vertical  velocities  along  the  £  axis  to  the  velocity  potential 
along  the  %  axis.  Thus,  we  require  that 

1 

v(£)  =  (ka)  ({,  £)j  I)  w  (I)  d|  (2.9) 

o 

for  0  5  £  - 1 }  since  $  (-£*  Oj  fc)  58  -  •  (t>  0}  £),  Eq.  (2.9)  is  auto¬ 
matically  satisfied  for  negative  {  ,  with  w(-£)  =  -  w(|). 

The  problem  now  is  to  find  the  kernel  function 

G(t,  1)  5  ®  (t,  0;  I)  (2.10) 

and  this  is  readily  accomplished  by  conformal  mapping. 


The  successive  transformations  shown  in  Fig.  2-5  map  the  region  of  Pig.  2-4 
into  the  entire  f  plane  J  the  resultant  transformation  is 


■IK  Hi)]’ 


(2.11) 


where 


7  «  ~r  .  Now  let  f^  and  fg  correspond  to  z  »  7  and  -7  , 


respectively.  The  complex  potential  of  a  sink  at  f^  and  a  source  at  f-,, 
each  of  strength  2,  is 

f  „  f 

F  «  -  J  lQg  (  fTT^)  +  6  (I)  (2.12) 

where  g  is  an  arbitrary  function  of  |. 
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Substituting  Eq,  (2.11)  into  Eq.  (2.12)  gives, in  the  z  plane,  the 
corresponding  potential 


Nov  choose  g  so  that  F(0)  =  0}  then 


Finally,  G(g,  |)  =  Re  F  (z)  is 


o(g,  i)  =  -  -  log  .(1~^)y(1-t'§)y,. 

"  (i+oy(i+T)y-  (i-s/a-Tf 


(2.13) 


Note  that  G  has  a  logarithmic  singularity  at  {  =  J  .  Note  also  that 


G(|,  I)  =  G(|,  |)  =  -  G(-|,  1) 

To  recapitulate,  then,  the  solutions  of  the  homogeneous  (singular)  integral 
equation 

v(l)  =  \\a)  j' G(g,  g)v(T)  d J  (2.14) 

o 

with  G(|,  |)  given  by  Eq.  (2.13),  provide  the  characteristic  free-surfaee 
shapes  v(g),  and  the  associated  eigenvalues  of  (Xa)  yield  the  natural 
frequencies. 
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The  numerical  solution  of  the  integral  equation  (2.14)  has  been  performed  for 
various  depths  on  the  basis  of  the  matrix  formulation  described  in  Appendix  A. 
The  results  are  presented  in  Section  4. 

Special  cases.  Several  special  cases  for  which  approximate  solutions  can  easily 
be  obtained  are  of  particular  interest.  For  the  half- full  canal  (e  =  0),  the 
kernel  becomes 


G(l,  I)  »  -  “  1°6  f|4)fl+il) 


A  one-term  Galerkin-type  solution  for  the  lowest  frequency  can  be  carried  out 
without  difficulty  for  this  case. 


Proa  Eq.  (2,14), 

1  11 

J  fw(5))2  <4  =  (Xa)  I)  w(I)  w(£)  d|dl  (2.15) 


Introducing  the  approximate  mode  shape  w  =  |  into  Eq.  (2.15)  and  carrying 
out  the  integrations  give,  for  e  =  0,  (a  =  l)  ,  the  approximate  fundamental 
eigenvalue 


h  "  Ts  71  -  1-367 

3J'8 


(2.16) 


J 

«  1.169  (|) 


Thi's  is  precisely  the  result  found  by  Rayleigh  in  an  entirely  different  way. 
The  result  is  necessarily  an  upper  bound  to  the  exact  frequency.  Rayleigh 


also  obtained  an  improved  approximation  a>  1.1644  (§) 

R 


e>V- 
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Also  of  special  interest  is  the  limiting  case  of  the  nearly-fuil  canal,  where 
0  i  and  7  -*  1.  The  integral  equation  becomes  simply 


w(£)  = 


(2.17) 


Denote  the  critical  values  of  (Xa)  by  c  (n  s  1,  2,  ,  .  .  )j  then,  for  e 
sufficiently  close  to  1  , 


c 


Thus,  the  frequencies  tend  to  infinity  as  the  full  condition  is  approached. 

This  result,  as  well  as  the  form  of  the  integral  equation  (2.17),  could  have 
been  anticipated  by  consideration  of  the  sloshing  of  liquid  in  a  half  plane 
with  a  rigid  boundary  having  an  aperture  of  length  (2aR).  (See  Fig.  2-6.) 

The  kernel  function  for  this  case  can  be  written  directly  in  terms  of  a  source- 
sink  pair,  without  the  need  for  conformal  transformations.  A  one-tera  Galerkin 
solution  of  Eq.  (2,17),  again  with  v  =  |  ,  gives  the  approximate  result 


so  that 


for 


e±  «  2*/3  =  2.094 


e  -*  1. 

■ 

Fig-  2-6  Limiting  Case  of  Nearly-Full  Container 


2.094 

JT7T 


(2.18) 


LOCKHEED  AIRCRAFT  CORPORATION 


2-10 


MISSILE  SYSTEMS  DIVISION 


faaMaunjwtKi  wbujhi  wnjBif i  WKW#*^ t» 


I24SD-5151 


For  the  nearly-empty  canal  (e  -*-l),  7  becomes  infin’.te.  The  asymptotic 
results  for  this  case  are  most  readily  found  by  abandoning  the  integral 
equation  and  exploiting  shallow-water  theory.  (See  Ref.  1,  p.  291.)  The 
differential  equation  for  the  velocity  potential  may  be  conveniently  deduced 
from  the  variational  principle  given  by  Lawrence,  Wang,  and  Reddy  (Ref.  J). 
The  variational  equation  is 


*[/H2«-£/  ♦*«]-<> 

L  V  F 


(2.19) 


where  V  denotes  the  volume  of  fluid  in  Fig.  2-1. 

As  shown  *n  Ref.  3>  Eqs.  (2.3)>  (2.4),  and  (2.5)  constitute  the  Euler  equation 
and  natural  boundary  conditions  of  Eq.  2.19. 


For  the  two-dimensional  canal  of  shallow  depth  variation  h(x),  it  is  assumed 
that  ♦  may  be  taken  independent  of  y.  Then  Eq.(2.19)  beccmes 

2  2 

t  V  I 

c 

■  IOXJ  dg  j  j 

F  F 


6 [l  /h(x)[i]  “-i  I'***]’ 


The  Euler  equation  is 


g 


(2.20) 


For  the  nearly  empty  circular  canal,  the  depth  has  approximately  the  parabo.lic 
variation 


hW  -  S  L 


(aS)S  -  x2J 


so  that  Eq. (2.20)  becomes 


s]*(4) 


❖  =  0 
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or,  in  terms  of  |  = 

I  —  0 

whence  p 

(r  •  1)  ♦"  +  2  {  ♦*  -  a  ♦  C  0  (2.21) 

This  is  essentially  Legendre's  differential  equation.  For  the  anti symmetrical 
modes,  $(0)  =  0;  in  addition,  the  condition  of  finiteness  it  {cl  is  imposed. 
The  solutions  of  Eq.  (2.21)  are  then  the  odd-order  Legendre  polynomials  of 
the  first  kind 


P2n-1  (*)  n  -  1,  2,  J,.  .  .(2.22) 

with  the  corresponding  eigenvalues 

X  =  2n2  -  n  n  .  1,  2,  3,.  .  .(2.23) 

n 


The  polynomials  (Eq. [2.22] )  represent  the  limiting  free- surface  shapes  of 
the  natural  inodes  as  the  depth  approaches  zero. 


The  first  eigenvalue  is  =  1  ,  and  so  the  fundamental  frequency 


is,  as  could  he  expected,  simply  the  frequency  of  a  pendulum  of  length  R  . 
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-.2*3  &maxu&  mx. 

*  ^  ’  -\ 

■'  The  gecigetrical  pararntters  of  Figs.  2-2  and  2-3  may  be  considered  to  apply  to 

V*^-\  *  -  *  v  - 

the  Spherical  tank.  A  rigorous  analysis  was  executed  for  only  the  nearly- 
fiOl  and  half-full  cases j  rigorous  resets  for  tin  nearly  empty  case,  given 
by  Lamb  (Ref.  1,  p.  291),  are  also  available. 

Introduce  the  cylindrical  coordinate  system  (r,  y,  0),  u^ere  x  =  r  cos  0, 
z  a  r  sin  C.  Hie  appropriate  nondim ens i onal  system  becomes  (p,  q,  S), 
where  p  «=  r/a R,  and,  as  before,  q  =  y/aR.  In  these  nondimensional  coordinates, 
the  governing  differential  equation  is: 


a2»  1 

bp2  p 


do 


b2* 

b^ 


L  *£± 
?  be2 


=  c 


(2.24) 


and.  the  freo-surface  boundary  condition  remains  that  given  in  Eo.  (2.7).  It 
is  evident  that  solutions  for  the  modes  may  be  written  the  form 


♦  *  t(p,  n)  cos  k  0 


(k  =  £,  2,  •  .  • ) 


(2*25) 


It  will  be  shown  later  that  only  the  modes  with  k  «  I  are  induced  by  lateral 
accelerations  in  the  5  direction.  Consequently,  attention  will  be  restricted 
to  these  cases. 


The  procedure  for  setting  up  an  integral  equation  analogous  to  Eq.  (2.14)  is 
as  follows.  Consider  the  volume  of  fluid  enclosed  by  the  wetted  surface  and  its 
reflection  about  the  originally  free  surface  q  =  0  .  (See  Fig.  2-4.)  Now, 
in  the  plane  q  =  0  (Fig.  2-7),  place  a  dis+rlbution  of  three-dimensional 
sinks  along  the  annulus  (p,  p  +  dp)  of  strengths  2f  (p)  cos  a  per  unit 
area.  Let.  f(p)  fl  (p,  q,  0}  p  )  p  dp  be  the  re  stilting  velocity  potential. 
Then,  since  the  vertical  velocity  just  below  the  plane  q  =  0  is,  because 
of  this  potential,  aero  everywhere  except  in  the  annulus,  where  it  is 
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f({5)  cos  9,  the  imposition  of  the  boundary  condition  Eq.  (2.7)  at 
9  =  0  requires  that 

1 

f(p)  =  (Xa) J  n  (p,  o,  o;  p)p  f(p)  dp 
o 

It  is  convenient  to  let  g(p)  =  Jp~  f(p)  j  then 

1 

g(p)  *  (\a)  /  H(p,p)  g  (p)  dp 
o 

where 

H(pip)  =/pp  0  (o,  0,  0;  p) 

The  problem  now  id  to  determine  H (p>  p)  . 


Fig.  2-7  Coordinates  in  Ere-  Surface  of 
Spherical  Tank 
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2.3.1  Nearly-Full  Tank 


For  the  nearly- full  case,  the  symmetrised  region  is  the  entire  space 
(Fig.  2-6),  and  the  potential  of  a  sink  in  the  plane  tj  =  0  is  unmodified 
by  the  presence  of  a  rigid  boundary  in  this  plane.  The  potential  of  a  sink 
of  strength  2  at  (p,  0,  «}  is 


1_ 

2* 


JTW 


_  2  —  2 
(p  cos  0  -  p  cos  a)  +  (p  sin  0  -  p  sin  a) 


and  so 


2* 


0  (p,  0,  Oj  p)  =  p— ■ 

o  t*3 


cos  a  d  a 

2""'“ - ~~Tl72 

J  -  2pp  cob  a j  ' 


Then, 


it 

/I  \  rr-  f  cos  a  d  Of 

*(p* p)  c  U !/*  J  r2~r~r 

v  Q  [p  +  p  -2  pp  cos  aj  • 


Clearly,  H  (p?  p)  =  H(p,  p).  Consider  p  <  p  ,and  let  m  = 


Then, 


H 


_  vZL.  f  _ cos  a  da 

*  J  jm2  -  2m  cos  a  +  lj  1^2 


(2.29) 


Letting  p  -  —  yields 


H  ■ 


rr 

2 


a  f  (g  °n.  s  - 1)  ap 

•j  [TT^Tfw 


(2.30) 


where 


2  JT 
q  =  bTTT 
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This  may  be  evaluated  in  terms  of  complete  elliptic  integrals  (Ref.  k,  p.  73) 
to. give 


(2.31) 


and 


2  4  2 
q.  sin 


d0 


are  the  complete  elliptic  integrals  of  the  first  and  second  kind,  respectively. 


The  result  (Eq.  [2.31]}  for  H  (p,  p)  may  be  put  into  a  much  more  convenient 
form  by  exploiting  the  Landen  transformations  (Ref.  5) 


K  <*■>=■  (!=)*<«> 

E  (■)  -  (s)  +  (^)k  (l) 


to  give,  as  a  final  result  ~ 

alp^)=  ff  HiHf)] 


(2.32) 
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for'  p  <  p  .  Since  H  (p,  p)  =  H  (p,  p),  the  kernel  is  canpletely  defined. 
This  kernel,  like  the  ones. for  the  circular  canal,  also  has  a  logarithmic 
singularity  at  p  »  p  . 

2.J.2  Half-Full  Tank 

For  the  half-full  case,  the  symmetrized  volume  is  a  complete  sphere,  and  so 
-■  the ‘kernel  function  can  be  evaluated  on  the  basis  of  the  known  Neumann  function 
(or  Dirlehlet  function  of  the  second  kind)  for  the  sphere.  Thus  (see  Bef.  6) 
the  function  of  the  position  Q  (Fig.  2-8)  given  by 

N(Q)  1  —  +  - 1 -  +  log  - aca — —r - - - — . - - — 

2*  )  PQ  (®)(QP*)  1  ~  (OQ)(OP)  cos  7  +  (OP)(QP*) 

satisfies  the  following  conditions  in  the  unit  sphere: 

(1)  N(Q)  is  harmonic. 

(2)  N(Q)  has  a  sink  of  strength  2  at  P  . 

(3)  The  nonaal  derivative  of  N(Q)  has  the  constant  value 
the  surface  of  the  sphere. 


J»* 


Fig.  2°8  Geometrical  Parameters  in 
Equatorial  Plane  of  Sphere 
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Nov  let  P  be  the  point  (p,  0,  0)  in  the  equatorial  plane  of  the  sphere 
(Fig.  2-7),  end  let  Q  be  the  point  (p,  0,  a).  Then 


and  therefore 


Note  that,  because  of  the  cosine  variation  of  local  sink  strength, at  the 
radius  p  ,  the  normal  flow  at  the  spherical  surface  vanishea^  as ‘--it-  rtwuld, 
despite  the  fact  that  the  Neumann  function  (Eq.  [2. 35])  for  a  single  slhk 

“  >■  7 

produces  a  constant  fLux  at  the  surface,  --**-* 
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which  gives  (Ref.  7,  p.  273) 


■s[ 


E(s)  -  (1  -  s2)  X  (s)1 


Then,  with  m  «=  2  (p  <  p),  evaluation  of  the  three  integrals  in  Eq..  (2.34)  gives: 
P 


H(p> p)  = [K(’>  -  E(”>] 


+  ~  K(b)  -  E(s) 


+  f= 

fT 


r[ 


E(s)  -  (1  -  s  )  K  (s) 


or,  finally. 


3/2  ; 

(pp)  x  (pp)! 


(2.36) 


for  p  <  p  .  Again,  the  additional  fact  that  H  (p,  p)  =  H  (p,  p)  '  r;_ 

(see  Eq..  [2.34])  completely  defines  H  (p,  p). 

As  in  the  case  of  the  circular  canal,  the  solutions  of  the  integral  Liquation  (2.27), 

with  the  kernels  Eq.  2.32  and  2*36,  were  performed  on  the  basis  of  • 

the  matrix  set-up  discussed  in  Appendix  A.  The  results  'are  given  in  Section  4.  . 

2.3.3  Nearly-Bnpty  Tank 

The  nearly- empty  tank  may  be  analyzed  in  a  fashion  analogous  to  the  treatment 
of  the  nearly- empty  circular  canal  by  approximating  the  wetted  surface  with 
the  paraboloid  of  revolution  having  the  depth  variation 


k  H  •  °] 
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and  using  shallow-water  theory.  The  problaa  of  oscillation  in  such  a 
container  has  been  solved  by  Lamb  (Ref.  1,  p.  291);  the  pertinent  parts  of 
the  solution  are  repeated  here  in  a  convenient  form.  By  assuming 
♦  «  f(p)  cos  6  (independent  of  q),  the  governing  differential  equation 
for  t  can  be  found  from  the  variational  principle  (Eq,  [2.19])  to  be 


p2(p2  -  1)  f*  +  p(3 P2-  1)  f  + 


£l  -  (2K  +  1)  p2Jf  - 


or 


|p(p2  -  1)  f  1  -  (2k  +  l)pjf  = 


(2.37) 

(2.38) 


Then,  appropriate  boundary  conditions  are  that  f  =  0  at  p  =  0  ,  and  f  is 

finite  at  p  =  1  •  Note,  then,  for  future  use,  that  Eq.  (2.38)  implies  that 

the  orthogonality  condition 

1 

fp  *n  (p)  tn  (p)  =  0  (2.39) 

o 

holds  for  eigenfunctions  f  and  f  associated  with  distinct  eigenvalues 

m  n 

X  and  X  . 
m  n 


The  solution  of  Eq.  (2.38)  for  the  antisyrametrical  modes  is 


where 


n 

-I  \  ^ 

J=1 


(2.40) 
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Twan  cnei 

U’war^v1, 


i,lj 

Ic 


l  ! 


C 


and  the  corresponding  eigenvalues  are 


\  a  2X1  -  1 

n 


n  s*  1,2, 


The  fi~st  three  eigenfunctions,  normalized  to  give  fn(l)  cl,  are: 


f2  «  5p5  -  2p 


f_  •  10p5  -  12p5  +  3p 


As  in  the  case  of  the  circular  canal,  the  lowest  eigenvalue  is  1, 

corresponding  to  the  frequency  of  a  pendulum  of  length  R  .  It  is  Interesting 
to  conpare  the  eigenvalues  for  the  nearly  empty  canal  and  spherical  tank; 
from  Eqs.  (2.41)  and  (2.23): 


^n^sphere  =  ^v)  canal.  +  (n  ’  -1-) 


LOCKHEED  AIRCRAFT  CORPORATION 


IMSD-5151 


Section  3 

SLOSHING  FORCES  AND  MOTION 


-  Siqspose  that  the  rigid  container  shown  in  Pig.  2-1  is  subjected  to  a  dis¬ 
placement  U(t)  in  the  x  direction;  tank  displacements  in  the  y  and 
%  direction,  as  veil  as  rotations  of  the  tank  about  all  axes,  are  assumed 
to  be  fully  constrained.  Define  a  displacement  potential  f  for  the  result¬ 
ing  fluid  displacement  relative  to  the  tank,  such  that  the  relative  displace¬ 
ments  are  given  by 

[$f  gf  d*  j 

dx  *  by  3  dz  J 

and  assume  ¥  in  the  form 

?  &n  ^  #n  V’  z ) 

-  where  t be  >  *3  are  the  potentials  associated  with  the  natural  modes  of 
n 

liquid  oscillation  in  a  stationary,  rigid  tank.  Lagrange's  equations  will 
be  the:  basis  for  establishing  differential  equations  relating  the  generalised 
coordinates  U(t)  and  a^(t)  to  the  external  force  X  on  the'  tank.  The 
results  vill  be  used  to  calculate  the  "sloshing"  forces  of  the  liquid  on  the 
container  vail. 


3-1 

LOCKHEH)  AIRCRAFT  CORPORATION  MISSILE  SYSTEMS  DIVISION 


The  potential  energy  during  motion  is 


-§/v[: 


ifv(l  ■.  & 


where  is  the  density. 
But  (see  Eq.  [2*5]) 


2 

co 

-~2-  ♦ 
g  n 
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4  P 


80  that 


1  p  /n  &  to  2  \i 


(3.1) 


Note,  however,  that  if  and  are  raodes  corresponding  to  the  distinct 

frequencies  w  and  w  , 
m  n 

f  7  *  *  7  ♦  d V  =  f  ♦  V2#  d V  -  f  ♦  7  *  •  iT  d  S 
Jo  n  Jan  Jo- 


(3.2) 
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Also ,  by  Eqs.  (2.5)  and  (2.4),  f  -  first  two  integrals  on  the  right-hand 
side  vanish,  and  so,  by  £4.  (2.5)>  the  right-hand  side  becomes 


*  dF 
n 


(3-3) 


But  then  the  left-hand  side  of  Eq.  (3*2)  is  also  given  by 


*  dF 
n 


(3.4) 


and  so,  for  m  *  o  ,  the  equality  of  Eqs.  (3*3)  and  (3^4)  implies  that 
*  m  a 


/  ♦  t  dF  -  0 

J  on 
¥ 


and  also 


•  V  ♦  dV  =  0 
n 


Hence,  the  potential  energy  (Eq.  [2.42])  becomes 


P.E. 


0> 


(3.5) 


(3-6) 
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or ,  letting  =  J  dF/ 


r.E.  " 


2g 


V  2  4 

)  a  a  v 

L  n  n  a 


(3.T) 


The  kinetic  energy  of  the  container- liquid  system  is 


2  PL 


K.E.  =  l/2  «  (U)  + 


TV  i  a*nl 

lL  n  TJ 


dV 


where  M  is  the  container  mass, 
c 


Using  Fqs.  (3.4)  and  (5.5),  together  with  the  energy  relation  for  a  single  mode 
(see  Eq.f2.l9]) 


gives 


2  P  P 

|7*n|  dV  "  ST/  *» 
F 


/  |VV  dV 

V 


dF 


K.E.  -  1/2  <MC  .  V  (U)2  +  ^  £  0^  «£  (i/  +  <3.8) 


where 


p 

«W  t! 


dV 


(3.9) 


and  is  the  total  mass  of  liquid. 

The  expression  for  0  may  he  transformed  as  follows: 

en  -M  f  ft)  +  e4  £  ft)  *4  f  ft)j  «  -/ 

v a*'  '  "  v 

ndS 


*?*♦„  av 

n 


S 
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Finally, 


2 

r 

3n  =  Tj  *  *n 


(3.10) 


With  the  Iagrangian  L-=(K«E.  -  ?.E.)  determined  hy  Eqs.  (3*6)  and  (3*7)> 
Lagrange * s  equations 


d_  dL  bL  Y 

atSu  -5u  =  x 


.  St  =  o 

at  aa„  “a 


yield  the  differential  equations 


a  +  J  a  = 
a  T  n  n 


-  (#.)  “ 


(3.H) 


<Mc  ♦  *L>  “  +  >t  = 


(3-32) 


The  sloshing  farce  Fg  of  the  liquid  acting  on  the  container  in  the  positive 
x  direction  is 


F0  =  M  0  -  X 
S  c 


Therefore,  from  Eq.  (3.12), 


PS  =  -  ,ml  u  "pl 


V  p  a' 

Z_,  &  n 


(3.13) 


These  results  vill  he  cast  in  convenient  forms  tor  the  circular  canal  and  the 
spherical  tank. 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILE  SYSTEMS  DIVISION 


lmsd-5151 


U4SD-5151 


(3.16) 


where 


ab  •/  t»n«)  ]2  *  i 


(3.17) 


and 


a. 

Bn  */  S  vn({)  4  ‘ 


(3.18) 


Here,  of  course,  all  forces  and  masses  are  understood  to  he  per  unit  length 
of  the  canal. 


It  is  Important  to  note  that  for  the  circular  canal  these  results  are  valid 
without  the  restriction  imposed  in  the  general  development  that  only  recti¬ 
linear  motion  of  the  tank  in  the  x  direction  was  permitted.  Rotation  of 
the  canal  about  its  center  produces  no  sloshing,  nor,  by  symnetry,  does 
vertical  tank  motion  produce  horizontal  sloshing  farce,  (in  fact,  it  can 
be  shown  that  vertical  motion  produces  no  sloshing  at  all  in  the  circular 
canal.)  Since  all  pressures  are  in  a  radial  direction,  the  resultant  horl- 
sontal  slosh  force  7g  always  passes  through  the  center  of  the  canal. 

3.2  SPHERICAL  TAM 


In  terms  of  the  sloshing  natural  modes  *  tJ1(r,y)  cos  k  0,  the  parameter 

VB 


0. 


^2  aS  2s 

**(*“  cos  Q)  (r#°)  eos  k  8  dr  &  0 


o  o 
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Thus,  it  follows  that  vanishes  for  all  k  ^  1,  Justifying  the  restric¬ 
tion  previously  made  to  consideration  of  only  the  case  k  *  1.  The  values  cf 
tn  at  the  free  surface  are  proportional  to  f^Cp),  where  fQ(p)  *  ~~  &n(p)t 

and  ^(p)  is  an  eigenfunction  of  the  integral  equation  (2.28).  Normalize 

fQ(p)  so  that  fQ(l)  *  1*  Then,  os  in  the  circular-canal  case,  the  slosh 

height  £  at  the  wall  is  a^c&^/g,  and 


pn  c  (aR)5  T~f f  p2  f(p)  eos2  6 


1  2s 


o  o 


2  1 
\5  -a 


*  (aR)^  f  P2  f(p)  *P 

o 


1  2s 


an  =  (aR)2  J  J  p[f  (p)l2  cos2  8  dp 


o  o 


=  s 


a- 

(sR)2  J  p  If  (p)l2  dp 


Hence,  from  Eqs.  (3.11)  -  (3«13), 

(M0  +  Mj  D  +  »pL  (aR)3£  Pnfs 


X 


Fs  -  -  K,  «  -  ,pL 


(*B)3  £  jy 


'a 
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Section  4 

NUMERICAL  RESULTS 


4.1  CIRCULAR  CANAL 


The  first  three  modes  and  frequencies  for  the  circular  canal  vere  found  for 
the  values  e  =  1,  0.5,  0.6,  0.4,  0.2,  0,  -0.2,  -0.4,  -0.6,  and  -0.8,  from 
the  20-by-20  matrix  formulation  described  in  Appendix  A.  The  values  of 
(^na)  and  XQ  (n  =  1,2, 3)  thus  obtained,  together  with  those  found  analyti¬ 
cally  for  e  =  -1,  are  given  in  Table  4-1. 


The  solid  curves  in  Fig.  4-1  show  the  variation  with  e  of  /\  (which  is 
proportional  to  frequency)  for  the  first  three  modes.  Worthy  of  note  is  the 
fact  that  the  variation  with  e  of  the  higher  frequencies  is  not  monotonic; 
the  minimum  frequency  of  a  given  higher  mode  appears  to  occur  slightly  below 
the  half- full  condition,  whereas  the  fundamental  mode  has  its  lowest  fre¬ 
quency  in  the  nearly- empty  state.  This  situation  is  at  least  partially  associ¬ 
ated  with  the  fact  that,  in  general,  as  depth  increases,  frequencies  of  fTee- 
surface  oscillation  tend  tc  shift  to  more  closely  packed  spectra.  Thus,  for 
example,  the  ratio  Xg/x^  decreases  monotonically  with  increasing  e  . 


M,  I  f! 


For  the  nearly- full  case  (e  =  l)  the  matrix  calculation  provided  the  result 
(4j»)  =  2.018  ,  as  compared  with  the  upper  bound  of  2.094  obtained  from  the 
one-term  Galerkin  solution.  The  corresponding  asymptotic  behavior 


(2.Q18)1/2 

(l-eS)1/1* 


for  e 


is  shown  as  -  the  dotted  curve  in  Fig.  4-1;  the  asymptotic  behaviors  for  Vx 
and  y  X^  are  similarly  indicated. 
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It  may  be  noted  that  the  upper  bound  =  1. 3^7  obtained  by  Rayleigh 
(or  by  the  one- term  Galerkin  solution  of  this  paper)  for  the  naif- full 
is  very  close  to  the  value  X^  =  1.3^0  found  from  the  matrix  cal¬ 
culation.  Rayleigh's  improved  upper  bound,  corresponding  to  X.-^  «  1.356, 
ia  actually  more  accurate,  but  there  is  little  doubt  that  the  results 
obtained  from  the  matrix  set-up  are  entirely  satisfactory  for  practical 
application. 

The  modal  parameters  and  defined  by  Eqs.  (3*16)  and  (3*17)  nre 
plotted  against  e  in  Fig.  4-2.  Except  for  the  nearly-empty  case  (e  =  -l), 
these  parameters  were  found  by  numerical  integration.  For  e  =  -1,  the 
result  (Eq.  [2.22]  )  for  the  nature!  mode  shapes  provides: 

A„  -  roc 

Bn  =  l/3  f or  n  =  1 
=  0  for  n  >  1 

th 

The  modal  parameter  B^  is  closely  related  to  the  amount  of  n  mode 

induced  by  lateral  acceleration.  The  low  values  of  B  for  n  >  1  indicate 

n 

that,  in  general,  the  higher  modes  would  not  have  a  major  influence  an 
sloshing  forces. 

The  mode  shapes  w q(£ )  are  shown  in  Fig.  4-5  for  e  =  1,  0,  and  -1.  A 
point  of  minor  theoretical  interest  in  connection  with  the  mode  shapes  is 
that  the  slope  at  §  =  1  can  be  shown  to  be  given  by 


Thus,  the  slope  is  negative  for  e  >  0  and  positive  for  e  <  0  j  for  e  =  1, 
the  slope  is  infinite.  The  nature  of  these  slopes  is  a  direct  consequence 
of  the  inclination  of  the  container  wall  at  the  free  surface. 
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4.2^  SPHERICAL  TANK  "  -  '  '  / 

The  values  of  (XQ a)  and  \Q  for.  the  nearly- full  (e  -  l)  and  half- full 

(e  *  0)  spherical  tank  were  found,  for  n  Ss  1,  2,  and  3j  free  the  matrix 
formulation  described  n  Appendix  A  and  are  given  in  Table  4.2.  Also 

presented  in  this  table  are  the  values  found  analytically  for  the  Dearly 

empty  tank  (e  =  -l) .  From  the  values  of  ( X^a)  for  e  =  1,  the  asymptotic 

behavior  ofvrn  near  e  =  1  is  shown  as  the  dot-dash  curves  in  Fig.  4-2,  and  I 

the  values  of  . /T '  for  e  =  0  and  e  =  1  are  shown  as  the  circles  in  this 

figure .  Toe  dashed  curves  represent  estimates  of  the  variation  of  in  ! 

the  intermediate  ranges,  based  on  the  presumably  analogous  trends  for  the 

circular  canal. 


The  modal  parameters  C  and  D  were  calculated  (by  numerical  int  gration 
“  n  n 

for  e  •«  1  and  e  =  0,  and  analytically  from  Eq.  2.42  for  e  =  -l)  and  are 
shon.i  in  Fig.  4-4-.  The  curves  connecting  the  calculated  points  in  this  figure 
were  estimated,  with  the  results  of  Fig.  4-1  for  the  canal  as  a  guide. 

Finally,  the  mode  shapes  i'n(p)  =  g^i'p)  for  the  three  calculated  cases 

\r  P 

are  given  in  Fio.  4-5. 
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TABLE  4-1 


EIGENVALUES  FOB  CIRCULAR  CANAL 


e 

a 

V 

\Qe. 

V 

H 

X2 

,  X3  7 

-1.0 

0.0 

0.0 

0.0 

0.0 

1.0 

6*0  .. 

15.0 

-0.8 

0.6  -- 

0.627 

3.23 

6.51 

1.045 

5.38 

10.85 

-0.6 

0.8 

0.879 

3.98 

7*30 

1.099 

4.97 

9.13 

-0.4 

0,917 

1.068 

4.34 

7.63 

1.165 

4.74- 

8.33 

-0.2 

0.980 

1.224- 

4.56 

7.82 

1.249 

4.65 

7-99 

0.0 

1.0 

1.3,60 

4.70 

7-96 

1.360 

4.70 

'7.96 

0.2 

0.980 

1.482 

4.81 

8.C6 

1.7X3 

4.91  - 

8.23 

0.4 

0.917 

1.596 

4.89 

8.15 

1.742 

5.34 

8.89 

0-6 

0.8 

1.706 

4.97 

8.22 

2.13 

6.22 

-  10,28 

0.3 

0.6 

1.822 

5.05 

8.30 

3.04 

8.42 

13.84 

1.0 

0.0 

2.018 

5.20 

8.44 

00 

00 

TABLE  4-2 

EIGENVALUES  FOR  SPHERICAL  TASK 


e 

w 

V 

v 

*1 

0.0 

0.0 

0.0 

0 

1.0 

1.565' 

5.34 

1 

0,0 

2.78 

5.99 

V 

>1 

X2 

*3  /, 

0.0 

1.0 

7.0 

17.0 

8-66 

1.565 

5.34 

8.66 

9.S5 

00 

<4 

ar-- 
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Section  5 

COMPARISON  WITH  EXPERIMENT 


Resonance  tests  (unpublished)  conducted  by  the  Lockheed  Missile  Systems 
Division  on  partially  filled  spherical  tanks  provided  the  experimental 
results  for  -/  ,  which  are  shown  as  crosses  in  Pig.  5-1*  25ie  estimated 

curve  for  given  in  Fig.  4-2  is  reproduced  in  Fig.  5-1  for  comparison 

She  agreement  is  generally  good,  and  the  agreement  of  experiment  with  the 
rigorously  calculated  frequency  for  the  half- full  case  is  excellent. 


Fig.  5-1  Experimental  Frequencies  for  Spherical  Tack 
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Section  6 
CONCLUSIONS 


The  results  for  the  circular  canal  are  sufficiently  complete  to  be  used  in 
the  dynamic  analysis  of  such  canals  filled  to  arbitrary  depth.  Some  degree 
of  caution  is  evidently  required,  however,  in  the  use  of  the  estimated 
results  for  the  spherical  tank  in  the  ranges  intermediate  to  the  nearly- 
empty,  half- full,  and  nearly- full  cases. 
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Appendix  A 

NUMERICAL  SOLUTION  OF  INTEGRAL 
EQUATIONS 


The  general  kernel  function  (Eq„  [iJ.13])  for  the  circular  canal  has  a 
logarithmic  singularity  at  -  g  >  G  .  For  |  near  |  ,  with  I  ^  1> 


where 


G{ it  I)  *  -  ~  log  |  i  -  I  |  +  p  (I) 


Ar-i 


p  (s)  =  -~iog  < c 

*  i (m)  7  -  d-5) 7 


(Arl) 


(A.2) 


K 


Near  {  =  1  , 


G(l,  5)  «  -  J  log  1 1  -  5  J  +  7  103  2 


(A*3) 


Note,  however,  that  G(0,  I)  =  G  ({,  0)  «  0  ,  and  recall  that  v(0)  *  0. 

1 

Then  the  integral  n  _  „  _ 

JG{ I,  5)  ’4  5)  as 

o 


u ay  be  approximated  by 


n 


^  CO  (1A)  J  G(£,  I)  al  -  w(l)  j  G(|,  I) 


(A.4) 
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vfcere  A  *  ~  .  Suppose  4  «  JA  ;  then 


^2 

j  0(4,  !/ 


dl  «  A  G(«5A,  1A) 


i  4  i  M 


0(5i  I)  dl  *  f  G(jA,  1) 


j  ^  N 


^1 


f  o  (6,  1)  dl  ««  Ap(jA)  +  f  -  -  log|t| 

“  A  JA  « 


~  A 
jA~2 


=  A^p(JA)  +  ~  (1  -  log  |)J  j  ^  N 
J  G(l,  I)  d5  -  f  [(  l  log  2)  +  2-  (1-  log  |)J 


Hence,  the  integral  equation  (2.14)  ia  approximated  by  the  matrix  equation 


AT  {v}  =  (  ^  )  fw} 


(A.  5) 


where'  {wj  is  the  column  vector  with  elements  v  (JA)  (j  «  1,  2,  .  .  .  »),  T  ie 


the  "integrating"  diagonal  matrix  r  1, 


j  and  A  is  the  syara^trical 
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c 

Si 


r 

* 


*»*» 

It 


square  matrix  with  the  general  element  A, ,  defined  by: 

i  k  i  k  N 
i  -  3  k  N 
i  «=  J N 

Note  that,  in  the  absence  of  a  logarithmic  singularity,  the  integration  scheme 
would  reduce  precisely  to  the  conventional  trapezoidal  rule. 

The  orthogonality  condition  for  distinct  eigenfunctions  w  and  w 

m  n 

1 

f  \  M  vn  ^  d*  “  0 
o 

is  approximated,  according  to  trapezoidal  rule,  by 

[wj  T  (o>n}  «  C  (A.6) 

where  [w^]  is  a  row  matrix  and  (w^)  is  a  column  matrix. 

The  orthogonality  condition  (Eq.  A.6}  is  satisfied  exactly  by  the  eigenvectors 
of  the  matrix  equation  (A.  5),  and,  for  this  reason,  the  present  integration 
scheme  ie  considered  to  be  a  reasonably  appropriate  one. 


Au  «  G(1A,  4A) 

«=  P  (JA)  +  j  (1  -  log  |) 

«  Jioe  2  +J  (1  -  iogf) 

7t  It  ^ 


f . 
1j 


L 


An  entirely  analogous  scheme  is  adaptable  to  the  integral  equation  (2.22)  for 
the  spherical  tank.  In  the  case  of  the  nearly-full  tank,  the  behavior  of 
H(p,  p)  near  p  =  p  is 


H  (p,  p)  «  -  ~  log  p  -  p  +  i  (log  p  +  log  8-2).  (A.T) 

71  71 


Including  the  case  p  =  1. 
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Eor  the  half- full  tank,  for  p  near  p  f  vith  p  !|  1, 


H  (p,  p)  m •-  j  log  j  p  -  p  |  +  -  [log  p  +  2p5  K(p2)  +  log  8  -  2] 


(A.8) 


X5n  the  other  hand,  H(l,  p)  near  p  *  i  is 


E  (1,  p)  «  -  |  log  (1  -  p)  +  j  (2  log  8  -  2) 


(A.9) 


Hence,  the  matrix  formulations  for  the  spherical  tank  may  "be  made  in  a 
fashion  similar  to  those  for  the  circular  canal J  the  eignevectors  are  then 
proportional  to  yp3F(p)  . 


Eigenvalues  and  eigenvectors  of  the  matrix  C  =  AT  (Eq.  [  A.  5  3  )  were  found  by 
using  a  Remington  Rand  1103AF  computer  on  the  basis  of  2Q-by-20  matrix 
approximation  (i.e.,  N  =  20).  The  highest  eigenvalue  of  C  (and  hence  the 
lowest  value  of  X  )  and  the  corresponding  eigenvector  were  fourd  by  matrix 
Iteration.  The  second  mode  was  found  frets  the  modified  ("swept ")  matrix 


C  *  C  •  C  (w1)  (w1l  T 


where  w  is  the  previously  determined  first  eigenvector,  normalized  according 
to 

[Vj]  T  {w1)  =  1 

After  determining  the  second  eigenvalue  and  eigenvector,  again  by  iteration, 

C*  was  similarly  swept  to  provide  a  matrix  C",  direct  iteration  of  which  gave 
the  third  eigenvalue  and  eigenvector. 
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